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Abstract 

We study a model quantum dot system in an external magnetic field by 
using both the spin-density-functional theory and the current-spin-density- 
functional theory. The theories are used with local approximations for the 
spin-density and the vorticity. The reliabilities of different parametrizations 
for the exchange-correlation functionals are tested by comparing the ensu- 
ing energetics with quantum Monte Carlo results. The limit where the vor- 
ticity dependence should be used in the exchange-correlation functionals is 
discussed. 
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I. INTRODUCTION 



The success of the spin-density-functional theory (SDFT)^ within the local-spin-density 
approximation (LSDA) to predict accurate results for the electronic structure of two- 
dimensional quantum dot systems depends on the exchange-correlation functional used. 
Until recently, functional based on the diffusion quantum Monte Carlo (DMC) calculations 
by Tanatar and Ceperley^ have been widely employed. The simulations were performed for 
the spin-compensated and spin-polarized two-dimensional electron gas (2DEG). Recently, 
Attaccahte et al.^'^ made fixed-node DMC calculations with improved accuracy for the 
2DEG including several partial spin-polarizations. They provided interpolation formulas 
which fulfil several exact results at the low and high electron density limits. 

In the current-spin-density- functional theory (CSDFT), used to describe electron systems 
in magnetic fields, the exchange-correlation energy depends also on the paramagnetic current 
density^. In the local approximation the current-dependence is converted to the dependence 
on the vorticity or on the Landau-level filling factor. The low fiUing factor (strong magnetic 
field) limit for the totally spin-polarized electron gas, is well known from the works by 
Levesque, Weis and MacDonald®, Fano and Ortolani^, and Price and Das Sarma^. The 
problem is how to interpolate the exchange-correlation energy between this limit and the high 
filling factor (zero-magnetic field) limit for a given spin-polarization. Several interpolation 
schemes have been suggested^^^^. 

The purpose of the present work is to study the reliabilities of the two zero- field exchange- 
correlation functionals and various exchange-correlation interpolation schemes with respect 
to the filling factor in the CSDFT calculations. We use a parabolically-confined quantum 
dot as the test system and compare our SDFT and CSDFT total energies with those ob- 
tained in variational quantum Monte Carlo (VMC) calculations^^'^^. In the zero-field case 
the parabolic confinement is lowered towards the limit of Wigner-crystallization and the 
difference between the total energies of the spin-polarized and spin-compensated solutions 
are monitored. For a finite confinement we calculate the ground state energy as a function 
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of the external magnetic field. Besides the CSDFT schemes, calculations are performed 
also with the SDFT, ie ignoring the current-dependence of the exchange-correlation energy. 
Thus, our calculations enlighten how important improvement the CSDFT is in comparison 
with the SDFT. For the quantitative comparisons we need numerically accurate results. We 
perform the SDFT and CSDFT calculations on two-dimensional point grids imposing no 
symmetry restrictions^^'^^. The convergence with respect to the grid size and other numeri- 
cal approximations is carefully tested. The VMC results are converged beyond the statistical 
noise. 

The outline of the present paper is as follows. In Sec. II we shortly describe the exchange- 
correlation functionals in zero magnetic field and give our results for a six-electron quantum 
dot as a function of the confinement. In Sec. Ill we discuss the interpolation as a function 
of the Landau level filling factor and show our results for the magnetic field dependence on 
the ground-state total energy. Sec. IV contains the conclusions. Effective atomic units are 
used in the formulas throughout the paper and for presenting the results they are converted 
by using the material parameters of GaAs, i.e. the dielectric constant e = 12.4, the effective 
mass m* — 0.067 and the gyromagnetic constant g* — —0.44. We choose the coordinate 
system so that x and y are in the plane of the dot and the z-axis is perpendicular to the 
plane. 

II. ZERO-FIELD AND THE LOW-DENSITY LIMIT 

Attaccalite et al. (AMGB)^'^ calculated the ground state of the 2DEG with the fixed- 
node DMC. Compared to the work by Tanatar and Ceperley (TC)^ there are a number 
of improvements in the numerical calculations. Backfiow correlations in many-body wave 
functions are included, and infinite size extrapolations are performed in the Monte Carlo 
data. An important new feature of the results by Attaccalite et al. is the appearence of a 
spin-polarized ground state before the Wigner crystallization. 

On the basis of their DMC calculations, AMGB^ proposed a new analytic representation 
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of the correlation energy. It takes into account several exact results in the low and high 
density limits. They give the exchange-correlation energy of the homogeneous electron gas 

as 

exc(r„ C) = e,(r., C) + {e'^''' - l)ef (^., C) + «o(r.) + a,{rs)C + «2(r,)C^ (1) 

Here = Xj-sfwn is the density parameter, calculated from the two-dimensional electron 
density n, and C, — {n^ — n])/n is the spin polarization determined by the spin-up and 
spin-down electron densities n-j- and nj^, respectively, exirs-, C) is the exchange energy, which 
can be written as 

ex(r„ C) = -2v^[(l + 0=^/' + (1 - C)'/']/37rr,. (2) 

In Eq. (1), ei^ contains the terms of the Taylor expansion of Cx with respect to ^ at ^ = 
which are beyond the fourth order in C, ct's are density dependent functions of the generalized 
Perdew-Wang form^^, and /3 = 1.3386. 

The above form (1) is based on DMC calculations for partial polarizations (0 < ^ < 1). 
TC^ made DMC calculations only for spin-compensated (C = 0) and spin-polarized (C = 1) 
systems. To calculate the correlation potential for finite polarizations the TC data is often 
used with the exchange-like interpolation^^ of Eq. (2). This leads to deviations from the 
exact results at high electron densities and too small spin-susceptibilities at low densities^. 

Our test system is a six-electron quantum dot confined by a parabolic external potential 

Kxt = ]^'rn*ulr'^, (3) 

where fiuiQ is the confinement strength. First we use a zero magnetic field and compare 
the results obtained with the two LSDA functionals to the VMC results. We study the 
weak-confinement limit at which the electron density in the dot is low and therefore the 
contribution of the exchange-correlation energy to the total energy is relatively large. Thus, 
the high-correlation effects in this particular system, i.e. the Wigner crystallization, spin 
density wave (SDW) formation, and spin-polarization, are assumed to be sensitive to the 
LSDA functional used. 
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Fig. 1 shows the energy difference between the Sz — 3 and Sz — spin states as a 
function of the confinement strength. The agreement with the VMC results is remarkably 
better for the AMGB functional than for the TC parametrization. In the AMGB results 
the spin-polarization is favoured at confinement strengths below hcuo — 0.23 meV. This is 
roughly in the middle between the TC and VMC results of 0.18 and 0.28 meV, respectively. 
As the confinement is made stronger, the difference between the parametrizations grows 
further enhancing the advantage of the new LSDA so that it smoothly follows the VMC 
curve. The origin of the difference is in the total energy of the polarized {Sz = 3) state. It 
is lowered when the AMGB functional is used instead of the TC one. This observation is in 
accord with the result of Gori-Giorgi et al^ that the improvement due to the new functional 
is directly proportional to the polarization and the electron density of the system. In the 
present case, however, there is a significant difference between these two approximations 
down to hcuo — 0.15 meV. This corresponds to quite a small electron density, i.e. rg — 
[N^/^Uloq)'''^^^ ~ 14 (Ref. 16). In our calculations, this is actually the smallest confinement 
strength for which we have strictly converged SDFT results. 

The transition point to the -S";;; = SDW state also seems to depend strongly on the 
applied type of LSDA. By using the TC functional we find the breaking of the spin symmetry 
at fiujQ ~ 0.45 meV, whereas with the new LSDA the transition occurs already at Tiujq ~ 0.8 
meV, which corresponds to = 4.5. This should be compared with the estimation by 
Koskinen et al. that in the case of closed shells, a SDW is found for > 5 (Ref. 16). Fig. 
1 shows also the results obtained by forcing the spin densities to be equal, i.e. preventing 
the SDW formation. It can be seen that the values of the SDW solutions are in a clearly 
better agreement with the VMC results than those of the symmetry-restricted solution. The 
relative amplitude of the SDW grows rapidly at low confinements below the transition point, 
and the electron density shows localization around six maxima, corresponding to the Wigner 
crystaUization. The behavior of the electron density is presented in more detail in Ref. 13. 
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III. EXTERNAL MAGNETIC FIELD 



Next we add an external magnetic field perpendicular to the dot (i.e. parallel to the 
z-axis) and fix the confinement strength to 5 meV. SDFT calculations are done with a 
minimal substitution of the external vector potential in the Schrodinger equation, and by 
using the zero- field exchange-correlation functionals for the 2DEG. In the CSDFT^ the 
exchange-correlation functionals depend on the electron currents in the system, and they 
are functionals of the spin densities and the vorticity. 



7(r) = V X 

n{r) 



(4) 



Above, n is the electron density and jp is the paramagnetic current density. Equally, we can 

use the Landau level filling factor given by 

v — 27771/7. (5) 

The data for the 2DEG in magnetic field is scarce. Fano and Ortolani^ used the Monte 
Carlo data by Morf et al}"^ and gave the v dependence of the 2DEG exchange-correlation 
energy for low v values from to about 0.8, corresponding to high magnetic fields. At the 
limit of the infinite magnetic field (z/ ^ 0) , the resulting curve agrees well with the results 
by Levesque, Weis, and MacDonald^. Price and Das Sarma^ obtained for the polarized 
2DEG the density dependence of the exchange-correlation energy at the v values of 1/7, 
1/5, 1/3, 1 and 2. In the low v region the fractional quantum Hall effect causes cusps in the 
exchange-correlation energy. Heinonen et al}^ took them into account in their modeling, 
but we have ignored them because in our test case of the six-electron quantum dot, the u 
values were above 0.9 in the magnetic fields up to 10 T. 

In the present calculations we have used expressions given by Koskinen et al}^ and 
Ferconi and Vignale^^ for the magnetic-field dependence of the exchange-correlation (See 
also Ref. 9). Koskinen et al. fitted their functionals in the high vorticity limit and used the 
following formula for the interpolation to the zero- field limit, 

e^,(n, C, v) = -0.782V2^e-^W + =°(n, C)(l - 6"^^), (6) 
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where /(i/) = l.di/ + 7u^. Koskinen et al. used the TC functional as the zero-field limit 
e^^^. We replace it by the AMGB functional (Eq. (1)). Ferconi and Vignale apphed a Pade 
approximant, fitting the low u limit of Levesque et al.^ to the zero magnetic field functional, 
i.e., 

elt''''{n, C, = (e™(n, u) + u^e^^'i^, C))/(l + (7) 
where e^^^^ is the interpolation formula for the infinite magnetic field limit^, 

gLWM ^ _o.782133V27m(l - 0.211i/°-^^ + O.OUu^'^). (8) 

As before, we insert the AMGB functional for e^c^^- The above formulae interpolate between 
the fully polarized electron gas values at high magnetic fields and the zero- field limit, which 
may have arbitrary polarization. Data for the intermediate polarizations at high magnetic 
fields would be desirable to test the interpolation further. 

Fig. 2(a) shows the quantum Monte Carlo data by Price and Das Sarma^ and Fano 
and Ortolani^ for the exchange- correlation energy of the polarized 2DEG at n = 0.138 a.u. 
The horizontal fine indicates the zero-field hmit of AMGB. Fig. 2(a) shows also values from 
exchange-correlation energy by Koskinen et al. (Eq. (6)). At low u values all the data agree 
well but when u increases the results by Price and Das Sarma (i/ = 1 and u = 2) approach 
much more slowly the zero- field hmit than the data by Fano and Ortolani. Therefore we 
propose a new functional which is an interpolation to the data by Price and Das Sarma. 
The resulting functional is 

e^^"(n, C, 1^) = (e™(n, i/) + i^^ef =°(n, C))/(l + 0.0061i/ - 0.03141/^ - 0.02011/^ + i^*), (9) 

and it is shown in Fig. 2(b). In the same figure we have also plotted e^^^^ of Ferconi and 
Vignale. One should note that at < 1 it approaches the zero-field limit much more slowly 
than the data by Fano and Ortolani in Fig. 2(a). The e^^^^ shows a maximum above the 
zero-field values, for which there is, however, no physical reason. 

The exchange-correlation potential V^c is obtained in the local approximation by calcu- 
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lating the functional derivative 

V^cA'^h ^) = d{ne^c)/dna. (10) 

In the CSDFT, the vector potential depends on the derivative with respect to the vorticity. 
The X and y components of the exchange-correlation vector potential Axe are 

A.e. = l/n#^^, (11) 



xc. = -yn-^^. (12) 



d d{ne^ c) 
dx d'j 

It should be noted that the contribution of the exchange-correlation energy to the total 
energy arises via the value of exc at the given density and vorticity and via its derivatives 
with respect to both density and vorticity. Therefore the total energies are sensitive to even 
slight variations in the exchange-correlation functional . 

The vorticity is a functional of the ratio between the paramagnetic current density and 
the electron density (Eq. (4)). Its values increase rapidly towards low density areas causing 
numerical instabilities. To circumvent them we have used the convoluted form of Koskinen 
et al}^ in the calculation of the partial derivatives de^c/d^ in Eqs. (11) and (12), 

^^''-'l=^e..inXn'). (13) 



d-f J V2^A3 

Above, the width A of the Gaussian function should be carefully adjusted to the vorticity 
of the system at the given magnetic field strength. Too small values for A may result in 
convergence problems. On the other hand using too large a value for A may cause inaccuracy 
in the results. In our calculations A ranges from 0.025 to 0.05. 

Fig. 3 shows the total ground state energy of a six-electron dot as a function of the 
magnetic field up to 10 T. Results obtained with the SDFT and CSDFT formalisms and 
different functionals are compared with the VMC data. Assuming again that the VMC 
results follow most faithfully the exact ones, the CSDFT formahsm is a clear improvement 
over the SDFT formalism: The CSDFT results obtained with the AMGB functional using 
the Fade approximant for the filling factor interpolation, are closer to the VMC values 
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than the SDFT values already at magnetic fields shghtly above 2 T. In this region Sz — 0. 
The improvement of the CSDFT formahsm is also clear above 5 T, where — S and the 
maximum density droplet (MDD) has been formed. The use of the AMGB functional instead 
of the TC functional improves the SDFT results in the regions of large spin polarization, 
i.e., around B — 4.5T where — 2, and in the MDD region. This conclusion matches again 
well with the analysis by Gori-Giorgi et a/.^. 

In the CSDFT, the ground state energy depends clearly on the chosen exchange- 
correlation functional. As shown in Fig. 3 the interpolation suggested by Koskinen et 
al. (Eq. 6) does not cause a big difference between the CSDFT and SDFT results. This is 
because e^c saturates quickly to the zero-field values at about v — 1 due to the exponential 
factors. The clear difference obtained by using the Fade approximant (Eq. (8)) with an 
increasing magnetic field results in a much better agreement with the VMC data. Second, 
the results obtained with the AMGB functional coincide well in the MDD region with the 
VMC ones, whereas those obtained with the TC functional (not shown in Fig. 3) are higher 
in energy resembling the situation in the SDFT calculations. Third, the new functional (Eq. 
(9)) based on the data by Price and Das Sarma seems to overcorrect the SDFT results. 
We conclude that despite the good fit to the Monte Carlo data for polarized 2DEG, the 
derivatives of Cxc seem to be poorly approximated by Eq. (9) as well as the spin compensated 
C = values for exc- 

In the MDD region of the VMC results the 3^ — 2 state constitutes the ground state 
in the CSDFT calculations for the exchange-correlations energies given in Eqs. (7) and 
(9). The corresponding energies are given in Fig. 3 as x-marks and dots, respectively. 
However, we believe that the Sz — 2 groundstate is an artifact due to the interpolation 
of the exchange-correlation energy between the high-field totally-polarized and the zero- 
field partially-polarized electron gas values. Due to lack of simulation data the polarization 
dependence enters the exchange-correlation functional only through the zero-field limits in 
Eqs. (6), (7), and (9). The high- field = 2 CSDFT ground states have the same value 
for the total angular momentum, Lz — 15, as for the Sz — S MDD-state, but the electron 
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at the innermost orbital has a flipped spin. Therefore, Hund's rule should favor the totally 
polarized 5"^ = 3 state. Further, we see that also around B = 4.5T, where both the CSDFT 
and VMC give the = 2 state, the CSDFT clearly underestimates the total energy. 

IV. CONCLUSIONS 

We have studied the electronic structure of a parabolically-conflned six-electron quan- 
tum dot using the spin- density-functional and the current-spin-density-functional theories. 
In particular the effects of weak conflnement and strong external magnetic fleld were con- 
sidered. Our main aim was to investigate, how reliably the different local approximations 
for the electron exchange and correlation follow variational quantum Monte Carlo results, 
which we consider as benchmarks. We have especially investigated the various interpolation 
schemes between the spin-compensated and totally polarized electron gas and the vorticity 
dependence of the exchange-correlation functionals. 

We flnd that the new LSDA functional by Attaccahte et al. gives in the zero magnetic 
field much better results for the total energy than the old form by Tanatar and Ceperley. 
According to our results, the same is also true in an external magnetic field for a six- 
electron quantum dot with a typical confinement of hcuo — 5 meV. For this system the effect 
of currents on the exchange-correlation energy becomes important in magnetic fields larger 
than 2 - 3 T. In our calculations for these high magnetic fields the local Landau level filling 
factor u is of the order of 0.9. The simulation data for the homogeneous electron gas in this 
regime is scarce, which hampers establishing of interpolation forms between the low and high 
magnetic field limits (large and small Landau level filling factors). The Fade interpolation 
by Ferconi and Vignale gives a very good agreement with the quantum Monte Carlo results 
in the case of spin-compensated and spin-polarized states of our system. The exponential 
interpolation form by Koskinen et al. underestimates clearly the effect of the magnetic field. 
Our results for the partially spin-polarised states in magnetic field are remarkably worse than 
those for the spin-compensated and spin-polarized states. Thus, our tests with the CSDFT 
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call for further simulations for the homogeneous electron gas in order to determine reliable 
interpolation forms for the vorticity and polarization dependence of the exchange-correlation 
functional. 
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FIGURES 




FIG. 1. Total energy difference between the 5^ = 3 and S'^ = spin states in a six-electron 
quantum dot. The results obtained with the SDFT using the LSDA functional by TC^ and AMGB^ 
are compared with the VMC data. The onsets of SDW solutions with decreasing confinement are 
marked by arrows. The dotted curve shows the symmetry-restricted results, corresponding to the 
LDA. 
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FIG. 2. Exchange-correlation energy per electron as a function of the Landau level filling factor 
for the spin-polarized 2DEG with n = 0.138 a.u. (a) Monte Carlo simulation data calculated by 
Price and Das Sarma^ (open circles), and the interpolation forms by Fano and Ortolani^ and 
Koskinen et al}^ (Eq. (6)). The zero-field limit by AMGB is shown by the dotted line, (b) Data 
of Price and Das Sarma^, the Pade approximant used by Ferconi and Vignale^^ (Eq. (7)) and the 
new fit (Eq. (9)) to the data of Price and Das Sarma e°^^. 
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FIG. 3. Ground state energy (minus 6 huh = 6 H^/lUq + where lOc = eB/m* is the 

cyclotron frequency) for the six-electron quantum dot in parabolic confinement with the strength 
hcuQ = 5 meV. (a) SDFT results using the AMGB (solid line) and TC (dash-dotted line) func- 
tionals. The VMC results are marked with the bold dashed line. The x-marks denote the CSDFT 
results obtained using the exchange-correlation functional by Koskinen et al. (Eq. (6)). CDW 
denotes charge density wave state which breaks the rotational symmetry, (b) CSDFT results using 
the exchange-correlation interpolations by Ferconi and Vignale (Eq. (7)) (bold solid line) and the 
new functional (Eq. (9)) (dotted line). The best SDFT results are repeated (solid line). The 
zero-field AMGB exchange-correlation functional has been used in these calculations. The x-marks 
and dots denote the CSDFT results for the Sz = 2 state obtained using the Ferconi and Vignale 
functional and the new functional, respectively. The z-components of the total spin and and the 
total angular momentum are given in parenthesis as (L^, Sz) and arrows mark the transition points. 
It should be noted that the (7, 1) state is missing in CSDFT. 
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